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Abstract 



We propose a model-independent metliod to determine the magnitude of the 
Cabibbo-Kobayashi-Maskawa matrix element \Vub\ from exclusive B and D 
decays. Combining information obtainable from B — > pii^i, B — > K*ui>, 
D — > piv^ and D — > K*ii'i, a determination of \Vub\ is possible, with an 
uncertainty from theory of around 10%. Theoretical uncertainties in the B — > 
K* ££ decay rate are discussed. 
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I. INTRODUCTION 



In the minimal standard model the couplings of the IV-bosons to the quarks are given in 
terms of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, Vij, which arises 
from diagonalizing the quark mass matrices. In the minimal standard model {i.e., one Higgs 
doublet) it is this matrix that is responsible for the CP nonconservation observed in kaon 
decay. A precise determination of the elements of the CKM matrix will play an important 
role in testing this picture for the origin of CP violation, and will constrain extensions of 
the standard model that make predictions for the form of the quark mass matrices. 

The present value of the 6 ^ m element of the CKM matrix, |Kfe| ^ (0.002 - 0.005) |l| 
arises from a comparison of the endpoint region of the electron spectrum in semileptonic B 
decay with phenomenological models. In recent years there has been a dramatic improvement 
in our understanding of the theory of inclusive semileptonic B decays pj-Q. It was shown that 
the electron energy spectrum, dT/dEe, can be predicted including nonperturbative strong 
interaction effects that are parameterized by the matrix elements of local operators between 
B meson states. For typical values of the electron energy Eg, the lowest dimension operators 
are the most important and the small nonperturbative strong interaction corrections are 
dominated by only two matrix elements, one of which is already determined by the measured 
B* — B mass splitting [0,^. However, for the semileptonic decay rate in the endpoint region, 
(m^ — m\))/2mB < E^ < (m^ — m^)/2mB, (where low mass hadronic final states are more 
important) the nonperturbative strong interaction corrections are large and an infinite set 
of nonperturbative matrix elements are needed. It has been shown that the same matrix 
elements determine the rate for B ^ Xg'^ m the region where the photon energy is near 
its maximal value 0. In principle, experimental information on i? ^ Xs7 can be used to 
predict the electron spectrum in the endpoint region of semileptonic B decay, leading to a 
model-independent determination of \Vub\- 

In this paper we propose a method for getting a precise model-independent value for 
\Vub\ using exclusive B and D decays. Our approach gives a value of \Vuh\ that (apart 
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from some very small factors) is valid in the limit of SU{3) flavor symmetry (on the u, 
d and s quarks) or in the limit of S'f/(4) heavy quark spin-flavor symmetry (on the c 
and h quarks). Consequently, the leading corrections are suppressed by factors of the small 
quantity {ms/rric — ms/rnh) — 0.1 or (ms/lGeV) ■ [as{rnc)/T^ — Q;s(mfe)/7r] ~ 0.01, and a 
determination of \Vub\ with a theoretical uncertainty of about 10% is possible. 

Semileptonic D K* i vi decay [I = e, /i) has been studied extensively and the form- 
factors which characterize the hadronic D K* matrix element of the weak current have 
been determined (with some assumptions concerning their shape) from the data. In this 
paper we denote the form-factors relevant for semileptonic transitions between a pseudoscalar 
meson containing a heavy quark, H, and a member of the lowest lying multiplet of vector 
mesons, by g^^~^^\ fi^^^) and a:^^'^\ where 

{V{p', e) I g 7^ Q \H{p)) = i g^''-''^ e^,,. e*^ (p + p'Y {p - p'Y , (la) 
{V{p\ e) I g7,75 Q \H{p)) = /(^"^^^ e; + (e* ■ p) (p + p'), 

+ aL^^^)(e*-p)(p-p')M, (lb) 



and £0123 _ _5qj^23 = 1. The sign of g depends on this convention for the Levi-Civita tensor. 
We view the form-factors g, f and a± as functions of the dimensionless variable y = v ■ v', 
where p = mu v, p' = my v', and = {p — p'Y — ^"h + "^y ~ 2m/^ my y. (Note that even 
though we are using the variable v ■ v', we are not treating the quarks in V as heavy.) The 
experimental values for the form-factors for D K* I vg are 

(z,^^*) 0.17 GeV-^ 

^^^ = - l + 0.63(,-l) ' 

g^---\y)^~ 0.51 GeV-^ (2c) 
^ l + 0.96(?/-l) ^ ' 

The form factor a_ is not measured because its contribution to the D K* i z/£ decay 
amplitude is proportional to the lepton mass. The minimal value of y is unity (corresponding 
to the zero recoil point where the K* is at rest in the D rest-frame) and the maximum value 



of y is (m|, + m\*) / {2m£) mx*) — 1.3 (corresponding to maximal K* recoil in the D rest- 
frame). Note that over the whole kinematic range 1 < ?/ < 1.3 / changes by less than 
20%. Therefore, in the following analysis of B decays we can extrapolate / with a small 
uncertainty to a somewhat larger region, which in what follows we take to be 1 < y < 1.5. 
The full kinematic region for 5— s>p£z/£isl<y<3.5. 



II. SEMILEPTONIC B ^ pii?i DECAY 



The differential decay rate for semileptonic B decay (neglecting the lepton mass), not 
summed over the lepton type i, is 



dy 

where r = nip/mB and S{ii) is the function 

s{y) = 



48 7r3 



fB~^P){y) (2 + ?/2 _ + 3^2 



(3) 



+ 4Re[af ^''^(y) f^^'^Xy)] m|r (y - r) {y^ - 1) 



+ 4 
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{y) m%r'{l + r' ~2yr){y' -I) 



f^-^Ky) [1 + d^^'-'Ky)] (2 + - Qyr + ?>r') 



(4) 



The function 3^^^''^ depends on the ratios of form-factors o'^'^'^^ / f^^^^^ and g^^^P^ / f^^~*p\ 
We can estimate S{y) using combinations of heavy quark symmetry and SU (3) flavor 
symmetry. Heavy quark symmetry implies the relations 

1 f m£,\^/'^ 



ar'"'(y) 



9 



{y) 



2 \mB 
KniB 



aJmr) 



-6/25 



(5a) 
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nit 
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nil, 



-6/25 
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{D-.K*) 



(5b) 
(5c) 



S'f/(3) symmetry implies that the B^ p+ form-factors are equal to the B K* form- 
factors and the B^ p° form-factors are equal to 1/V2 times the B K* form-factors. 




1 1.1 1.2 1.3 1.4 1.5 

y 

FIG. 1. The function S{y) defined in eq. (Q) as a function of the kinematic variable y = v ■ v'. 
The dotted vertical line corresponds to the kinematic limit for D K* lu^. 

In the limit where the the heavy quark Q has large mass, the matrix elements in eqs. (|l]) 
depend on uiq only through a factor of ^vtlh associated with the normalization of the heavy 



meson states. Consequently, for large mc, (a 



)/(« 



a 



order Aqcd/'^c, so we can set a 



{D^K*) 



-a 



{D^K") 
+ 



in eq. yielding 



\niB/ 



-6/25 



a 



(D^K-) 



) is of 



(6) 



Using eqs. (|5a|) , (|5^), (H) and SU{?>) to get the — > iug form- factors (in the region 
1 < y < 1.5) from those for D —>■ K*ii'£, given in eq. (@), gives S{y) plotted in Fig. 1. 
We use asinih) = 0.22 and a,(m,) = 0.39. In Fig. 2 we plot S^^-'P^y) and 6^^^^'\y) as 
a function of y. The latter function (which will be used later in this paper) is denoted by 
the dashed curve. Perhaps the largest uncertainty in our analysis for 6 comes from setting 



a 



{D-.K') 



—a 



If a 



-A a 



then eq. (H) gets multiplied on its right 



hand side by the factor (1 + m£)/ms)/2 + A(l — m£)/mB)/2. In Fig. 3 we plot S^^~*p^ and 
^(B^K*) for A = and 2. 

Note that 6 is fairly small, indicating that a~^^^'^ and g^^~*P^ make a small contribution 
to S{y) (in the region 1 < y < 1.5), so even significant corrections to eq. (^ will not have 



0.08 
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FIG. 2. The function S{y) as a function of the kinematic variable y = v ■ v'. The sohd curve is 

d(^^p\y), the dashed curve is d^^^^'^y). 
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a large impact on S{y). We can use our prediction for S{y) to determine \Vub\ from the 
B pii'e semileptonic decay rate in the region 1 < y < 1.5. Our predicted S{y), in 
Fig. 1, gives a branching ratio of 5.2 |V^bp for p'^ii^e in the region 1 < y < 1.5 

(corresponding to 16GeV^ < < q^^-^ = 20GeV^, which imphes Ei > 1.6 GeV in the 
B rest-frame). While such a model- independent determination of \Vub\ may eventually be 
superior to a determination from a comparison of the endpoint of the electron spectrum with 
phenomenological models there will be a sizable theoretical uncertainty associated with 
\Vub\ determined in this way from order SU{3) violation and order l/m^fe corrections 
to relations (|^) and (|]). What is needed to get a value for \Vub\ with smaller theoretical 
uncertainties is an improved method for determining \f^^'^^^\'^ (1 + 5^^^'^'^). 

Our method for determining a precise value for \ Vub\ is based on the observation that the 
"Grinstein-type double ratio" (fiB^p) / fiB^K^yi^fin^p) / f(D^K')^ is equal to unity in 
three separate limits of QCD (isospin violation is neglected here): (i) the limit of SU{3) 
flavor symmetry, 0, where the strange quark mass is treated as small compared with a 

typical hadronic scale; (u) the limit of SU{4) heavy quark spin-flavor symmetry, nib^c oo, 
where the bottom and charm quark masses are treated as large compared with a typical 
hadronic scale; {iii) the limit = m^, where the bottom and the charm quarks are related 
by an SU (2) flavor symmetry. Consequently, 



/ 



f 



f{D~.K*) 



o 



'm^_rris_ rUs a.;(mc) - ^^(^fe) ' 
nic mb' 1 GeV tt 



(7) 



We propose to extract a precise value for (1 + S^^^'^^) using 



/ 



(B^p) 



(1 + 5(S-P))= 



fiD- 



>P) 



(8) 



Multiplying by the ratio of D decay form-factors above cancels out SU{3) violation not 
suppressed by factors of the heavy quark mass in the most important part of the B — »• 
piue differential decay rate, i.e., the factor of \f^^^^^\'^, leaving an uncertainty from SU{3) 
violation only in 6. Since as we have argued, \6\ is likely to be less than 0.15, the effects of 
SU (3) violation in it can safely be neglected. The plots in Figs. 2 and 3 show the kinematic 



sources of SU{3) violation in S arising from the fact that the p and K* masses are not equal. 
There are also contributions from SU{3) violation in the ratios of the form-factors a+// and 

g/f- 

In principle, the form factor f^^^P^ can be obtained from experimental information on 
the Cabibbo suppressed decay D plu^. However, at the present time, the small branching 
ratio Dll Br(D+ — > p°/iz/^) = (2.0l;i;|) ^ 10~^ has made extraction of the form factor /(^"^^) 
too difficult. It may be possible in future fixed target experiments or at a tau-charm factory 
to determine f^^^f\ Assuming this can be done, the factor \ f^^~^^*)\^ (1 + ^{b^k*)^^ -g ^^^q 
remaining ingredient needed for a determination of (1 + 5'-^^^^) via eq. (H). 



III. RARE B DECAYS 



One avenue to find the factor \f^^~*^*)\'^ (\ + ^(^^^*)) ^ges the exclusive rare decays 
B K* a OT B ^ K* V u, which may eventually be studied at hadron colliders, or at B 



factories. The effective Hamiltonian for these decays is |Tl|-P 



4Gf 



V2 



(9) 



where p is the subtraction point (hereafter we set p = rrij, and do not explicitly display the 
subtraction point dependence), and the operators Oj are: 



Oi = 


{sLa 7m bLa) 


[clp 7^ Clp) , 




(10a) 


02 = 


{sLa 7m his) 


Cl/3 7^ CLa) , 




(10b) 


03 = 


{S Lai fib La) 


{ul(3 1^ Ulp) + . . 


• + {hp Y bLp) 


(10c) 


04 = 


{sLa Ifi bLp) 


(uLfSl^ULa) + ■ ■ 


■ + {bLf3 Y bLa) 


(lOd) 


05 = 


is La 7m bLa) 


{urp YURfi) + .. 


■ + (bRpl^ bRp) 


(lOe) 


06 = 


{sLa Ifi bLfs) 


{uRfsY URa) + ■■ 


■ + ibRpY bRa) 


(lOf) 


07 = 


(e/167r^) rrib 


{sLa,,bR)F''\ 




(lOg) 


08 = 


(fif/ievr^) nib 


{sLcr^u bR) 




(lOh) 




FIG. 4. Feynman diagrams whose contribution to exclusive rates is neither included in the 
form- factors, nor in the effective Wilson coefficient Cg. The black square represents one of the 
four-quark operators Oi — Oq. 




FIG. 5. Feynman diagrams whose contribution to exclusive rates is part of the nonperturbative 
matrix element of Cg Og. 

Og = (eVlGTT^) {SL^M {ir^) , (lOi) 
do = (eVlGTT^) (sljM ihV) , (lOj) 
On = (eVlGTr^sin^^iy) {sl^M [^7'^(1 - 75) i^] ■ (10k) 

For B K* ££ we need the matrix elements of Oi — Oq and Os at order and to all 
orders in the strong interactions, and the matrix elements of O7, Og and Oio to all orders 
in the strong interactions. Among the contributions to the B — > K* il matrix element of 
Oi — Oq are the Feynman diagrams in Fig. 4, where a soft gluon (with momentum of order 
k -C ^/^) connects to the qq loop. We are interested in the kinematic region 1 < y < 1.5 
which corresponds to a i£ pair with large invariant mass squared q^ between 14.5 GeV^ 
and 19GeV^. In this kinematic region we have found by explicit computation that the 
contribution of the Feynman diagrams in Fig. 4 are suppressed by at least a factor of kj \/q^ 
compared, for example, to the contributions of the diagrams in Fig. 5. In the region of large 
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(compared with the QCD scale and the mass of the quark q) the qq pair must "quickly" 
convert into the (color singlet) ^l pair and hence the coupling of soft, long wavelength gluons 
to the qq pair is suppressed at all orders in QCD perturbation theory. Similar remarks hold 
for the matrix elements of O^. This "factorization conjecture" implies that for B K* ii 
at large we can take the matrix elements of Oi — Oq and Og into account by adjusting the 
coefficients of O7 and Og by a calculable short distance correction. In the next-to-leading 



logarithmic approximation Cg is replaced by an effective Cg{y) coupling [|T3 



Here 



Cg{y) = Cg + h{z, y) (3Ci + C2 + SCs + C4 + 3C5 + C,)-\ /i(0, y) (C3 + 3C4) 

-\KI, y) (4C3 + 4C4 + 3C5 + Ce) + ^ (3C3 + C, + 3C5 + Ce) . (11) 



h{z,y) = -hnz + ^ + ^x (12) 



--{2 + x)^\l-x\{ 
9 




— it: 



^ ; for X = iml/q'^ < 1 



2 arctan(l/-\/a; — 1) ; for a; = Aml/q'^ > 1 , 
with h{0,y) = 8/27— (4/9) [ln{q'^/ml) — iir], and z = rric/mh, r = mK*/TnB- On the right 
hand side of eq. (0) q^ = m\ + — 2mBmK* y should be understood. Fig. 5 is now 
part of the nonperturbative matrix element of CgOg. Note that eq. ([TTl) differs from Ref. 
T3| , since the one gluon correction to the matrix element of Og is viewed as a contribution 



to the form-factors in our case. 

Using mt = 175 GeV, rUb = 4.8 GeV, = 1.4 GeV, as{mw) = 0.12, ^^(m;,) = 0.22 and 
sin^ 9]v = 0.23, the numerical values of the Wilson coefficients in the leading logarithmic 
approximation are Ci = -0.26, C2 = 1.11, C3 = 0.01, C4 = -0.03, C5 = 0.008, Cq = -0.03, 
C7 = —0.32. The operator Og does not contribute at the order we are working. Cg, Cig and 
Cii depend more sensitively on mt (quadratically for mt ^ my/)- In Table I we give their 
values for mt = 165 GeV, mt = 175 GeV and mt = 185 GeV. 

In eq. (|TT|) the second term on the right hand side, proportional to h{z,y) comes from 
charm quark loops. Since is close to 4m^, one is not in a kinematic region where the 
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rrit = 165 GeV 


mt = 175 GeV 


mt = 185 GeV 




4.17 


4.26 


4.34 


Clo 


-4.21 


-4.62 


-5.04 




1.40 


1.48 


1.57 



TABLE I. Coefficients of the Og — On operators at the scale nib for different values of the 
top quark mass. Cio is calculated in the leading logarithmic approximation, while Cg and Cn are 
calculated to next-to-leading order accuracy]^ 

perturbative QCD calculation of the cc loop (or factorization) can be trusted. Threshold 
effects, which spoil local duality, may be important. (In the kinematic region near = 
the charm quarks in the loop are far off-shell and eq. (O) should be valid. However, in this 



region we cannot justify using eq. (|Tl|) for the light quark loops.) Later we examine the 
sensitivity of the B —>■ K* 1 1 rate in the kinematic region of interest to cc threshold effects. 
For slightly lower values of (or equivalently for larger values of y) than we consider, such 
effects are very important. The rates for B K* J /if) — > K* il and for B ^ K* ip' ^ K* 
are much greater than what eq. (0) would imply. The latter process occurs with the ip' on 
mass-shell at ?/ = 1.6. 

The hadronic matrix element of O7 is expressed in terms of new hadronic form-factors, 
g± and h, defined by 

{V{p\ e) I qa,, Q \H{p)) = gf^''^ e,,,^ e*' {p + p' f + /"""^''^ e,.A. e*^ {p - p'Y 

+ e,,^„ {p + p'Y {p - p'Y (e* ■ p) , (13a) 

{V{p', e) I g Q \H{pY = I gi"^""^ K {p + p'), - e; (p + p'Y] 

+ i g^-^^^ [e* {p -p')^- e* {p - p')^] 

+ I [{p + p'Y (P - p'), -{p + p'), (P - p').] (e* ■ P) . (13b) 



*For Cg in the next-to-leading logarithmic approximation terms of order Og are sub dominant, 
since the leading contribution to Cg is order ln(m^/m^) ~ l/ag. 
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The second relation is obtained from the first one using a^'^ = | et^'^^'^ uxc, 75. The differential 



decay rate for B K* I i (not summed over the lepton type £) is 



dr(5 ^ ir*. 



a 



s 2 
vfT^ r 



\C,{y)\'S'{y) + \C^o\'S{y) 



(14) 



dy 2Att^ KAtt 

where S{y) is given by the expression in eq. (^), with the form-factors replaced by those 
appropriate for B — > K*, and r = rnK'/friB- S'{y) is obtained from S{y) via the replacements 



f 



f 



g(^^-*^'\ml-m'j,,)+g'_^-'^ >mUl + r'-2yr) A{y) , (15a) 



a 



(B^K*) . 

a\ + 



/i(^-*^')m|(l + r^-2|/r)-(7V 



(B^K*) 



My) 



(15b) 

(15c) 



where A{y) = 2m}jCT /[m\{l + — 2yr) Cgiy)]. Since C7 is small compared to Cg, it is 
convenient to rewrite the differential decay rate as 



dr(5 ^ K* 
dy 



247r3 



AttJ 



\C,iy)\^ + \C,o\ 



(16) 



X 



/(^-^*)(y) [1 + 5(^^^*)(Z/)] \^y' -l{2 + y'-6yr + 3r') [1 + A 



where A contains the dependence of the differential decay rate on C7. 

Unitarity of the CKM matrix implies that jV^t^Vtbl — |V^1K:6| (with no more than 3% 
uncertainty), so that once A(?/) is known, a value of (1 -|- ^i^^^*)^ can be deter- 

mined from experimental data on 5 ^ K* ii. To find A(y) we use the relations between 
the tensor and (axial-)vector form-factors derived for large mi, in Ref. 
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(B-^K*) 



+ 9 



f{B^K*) I 9 JB-^K*) 

*K*) r +^9 



TUB niK* y 



{B-*K*) (B^K") 

9+ -9- 

f^iB^K*) 



rriB 

-2m^^(^-^*), 



2mB 



(17a) 
(17b) 
(17c) 



Recent lattice QCD simulations indicate that these relations hold within 20% accuracy at 
the scale of the B mass |jl5[ . In the limit where mi, is treated as heavy, a^^^ ■* + a"^^^ ^ 
is small compared with a[^^^ ■* — a^^^^ \ so eq. ( |17cD can be simplified to 



^We correct some obvious factor-of-two errors in Q. 
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Note that a similar simplification for g^^^ ■* + g^^^ ^ is not useful, because in eq. (|15aj) 
+ g- IS enlianced by compared to (7+ g- 
Using eqs. (^), ([T7E| ) and (|T8D, A{y) is expressed in terms of Cy, 

Cg, Cio, g(B^K') f f(B^K') ^{B^K') ^ j(^B-.K*) _ Using eqs. d) and (§ to relate ratios of 

B K* form-factors to ratios of D K* form-factors, we find that in the kinematic region 
\ < y < 1.5, A(?/) changes almost linearly from A(l) ~ —0.14 to A(1.5) ~ —0.18. The 



value of A at zero recoil (using ~ m^) does not depend on the ratios of form-factors [16 

'4 Re [Q (79(1)] 4|C7p" 



A(l) = i 

|Co(l)P + |C 



10 



|2 



(19) 



1 — r (1 — r)^ 

Even though there are l/nic corrections to eqs. (^) and (^), they do not affect A(l). Fur- 
thermore, A is small compared with unity and has a modest y-dependence. Consequently, 
l/mc corrections to the y-dependence of A, and l/m^ corrections to A(l) can only have 
a very small impact on a value of (1 + ^'yB^K*)^^ extracted from the B — > K* ii 

differential decay rate using eq. (|T6|). 

Using the measured values of the D K* i vi form factors and the heavy quark symme- 
try relations in eqs. (§) and (g) to get |/(s^^*)|2 (1 + ^(B-i^*))^ together with \Vch\ = 0.04, 
Tb = 1.5 ps and a{m\Y) = 1/129, we find that eq. (|TB|) gives a branching ratio of 2.9 x 10^^ 



for B — s> K* i£ in the kinematic region 1 < y < 1.5. 

The largest theoretical uncertainties in using B ^ K*ii for extracting |/(-B-^*)|2 + 
come from order as corrections to the coefficients of the operators O9 and Oio and 
our treatment of the B K*ii matrix element of the four-quark operators. It is h{z,y) 
that takes into account the cc loop contributions to the matrix elements of the four-quark 
operators. 

A comparison with a phenomenological resonance saturation model [l7] gives an indica- 



tion of the uncertainties in the prediction for B —>■ K* £i that arise from the fact that the 
kinematic region we focus on is not far from DD threshold. In this regard we note that 
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M^(„) [GeV] 


r^(n) [GeV] 


Br(?/;(") - 






3.097 


8.8 


10-5 


6.0 


10-2 


^(2) 


3.686 


2.8 


10-^ 


8.4 


10-3 


^(3) 


3.77 


2.4 


10-2 


1.1 


10-5 


^(4) 


4.04 


5.2 


10-2 


1.4 


10-5 


^(5) 


4.16 


7.8 


10-2 


1.0 


10-5 


^(6) 


4.42 


4.3 


10-2 


1.1 


10-5 



TABLE II. Mass, width and leptonic branching ratio of the 1 cc resonances 

using factorization to estimate the B K* ^p^"^ K* £ I matrix elements of the four-quark 
operators {^p^'^^ is the ra'th 1 cc resonance) we find that in a resonance saturation model 
h{z, y) in the second term of eq. ( plD gets replaced byf] 



h(z v)^-K—^^ 



(20) 



where r^(„) and M^(„) are the width and mass of the n'th 1 cc resonance. Experimental 



values for these quantities and the branching ratios to i£ are given in Table II. In eq. (^OD k = 
2.3 e**^"" is the factor that the B — » J/ipK* amplitude, calculated using naive factorization, 
must be multiplied by to get the measured B — ^ J/ipK* rate. Since the magnitude of 
K is large, we do not assume that eq. (|20| ) has the same phase {i.e., = 0) as naive 
factorization would imply. Replacing h{z, y) in eq. ( |TTD by the expression in eq. (pOD results 
in an effective coefficient of Og that we call Cg. A measure of the deviation of this model 
for the cc resonance region from the expression in eq. (|ll|) is given by d{y) defined by 



+ \C,,\'' = {\C,{y)\^ + \C,,\^)[l + d{y)]. 



(21) 



•f-For not near the resonances, there are uncertainties associated with the dependence. In 
eq. (20) factors of not associated with the resonance propagator are set equal to the square of 
the resonance mass. 
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FIG. 6. The function d{y) defined in eq. (21) as a function of the kinematic variable y = v ■ v'. 
The sohd, dash-dotted and dashed curves correspond respectively to ip^ = 0, 7r/2 and vr. 

In Fig. 6 we plot d{y) for 1 < y < 1.5. Note that part of h{z,y) is associated with cc 
pairs at large virtuality, and so is reliably reproduced by QCD perturbation theory. In fact 
h{z,y) is scheme dependent, and so d{y) is only a very crude measure of the uncertainties 
that arise from being near the cc threshold. The solid, dash-dotted and dashed curves 
in Fig. 6 correspond respectively to ip^ = 0, 7i/2 and tt. This analysis suggests that the 
uncertainty associated with the charm threshold region has on average about a 20% effect 
on the B K* il rate for 1 < y < 1.5. 

The uncertainties, involving the DD threshold region and the order as contributions to 
Cg and Cio, can be avoided if the decay B K* vv can be studied experimentally. While 
this will be difficult, the large missing energy carried by the neutrinos in the kinematic region 
we are interested in may help . The differential decay rate for B K* v v (summed 
over the neutrino flavors) is 

dr(5 ^K*vv) 



dy 



16 7r3 V27rsin2^ 



a 



16 TT'^ 



a 



27r sin^ 9 



w 



m%r'\Cr,\'S{y) 



m| |Cii| 



X 



/(^-^^*)(l/) ' [1 + 5^''-^'\y)] {2 + y^- Qyr + Sr^) . (22) 
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The coefficient Cu depends on the top quark mass (see Table I). Once the top quark 
mass is known more accurately, the B —>■ K* v v differential decay rate provides a way 
to get (1 + (jC^^-f^*)) that, from a theoretical perspective, is very clean. Recall 

that the function is the analog of S^^^^^ that occurred in B piui semileptonic 

decay, but it depends on ratios of i? ^ K* form-factors that occur, instead oi B ^ p form- 
factors. It is plotted in Fig. 2 with the dashed curve using eqs. (|^) and (H) to deduce the 
ratios of form-factors a^~*^ and g^^~*^''^ / f'-^^^*^ from the D K*Ivi form- 

factors. S^^^^*^ (y) is fairly small, and so even though there is SU (3) violation in the relation 
between S^^^^*') and 5^^^^\ this does not introduce a large uncertainty in our prediction 
for + using eq. (|). Using eqs. (D and (D to get |/(^-^*)|2 (l + ^(5-^*)) 

from the measured values of the D — * K* I vg form-factors, we find that eq. (0) implies a 
branching ratio of 1.9 x 10^^ for B K* vv in the kinematic region 1 < y < 1.5. 

The difference in the factor \/y'^ — 1 (2-|-?/^ — 6?/r-|-3r^) for r = mp/niB and r = itik* /niB 
divided by their sum is less than 3% for 1 < y < 1.5. Therefore, it is a good approximation 
to rewrite eq. (^, using eqs. (^, (^) and (^), as 



dy S\Vt*M^ \ a\C^^\ J ml^ dy 



(23) 



If S'f/(3) violation in the y-dependence of the ratio of D decay form-factors in eq. (^) is 
small then we can also compare integrated B decay rates to get a precise value for \Vub\- 
Assuming that the shape of the form-factors / are well approximated by simple pole forms 
and taking the pole mass for fi^^^*) to be 2.5 GeV (corresponding to the D** mass) and the 
pole mass for f^^~*p^ to be 2.4 GeV (corresponding to the D** mass), we find that the ratio of 
D decay form-factors squared in eq. (^) varies by less than 0.5% over the range 1 < y < 1.5. 
It may be possible to get some model-independent information on the y-dependence of the 
ratio fiD~*p)/fiD'-K') using the methods of Ref. ^ . 

The D semileptonic decay rate is almost completely saturated by the K and K* hadronic 
final states. The heavy quark symmetry relations in eqs. (^ and (^ do not imply that 
the rare decay mode B — > z/ z/ (and also B —>■ X^ii when the effects of the four-quark 
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operators are neglected) is also saturated by these states in the kinematic region that overlaps 
with the D decay. For some of the D decay phase-space is small compared with m|), 
while the scaling relations in eqs. and (|^) hold for c and h quark masses treated as large 
with y held fixed. 

IV. CONCLUDING REMARKS 

In this paper we have explored the use of exclusive B and D decays to ob- 
tain a model- independent value of \Vuh\ with small theoretical uncertainties. Our 
method is based on the fact that the Grinstein-type double ratio of form-factors 
^j^(B^p)^j^(B-.Kny(^f{D^p)/f{D^K*)-^ is equal to unity in the SU{3) limit, and in the limit 
of heavy quark symmetry. A determination of \Vub\ with an uncertainty from theory that 
is less than 10% is possible using information obtainable from the decay modes B —>■ p i 
B ^ K*ui^, D ^ pi i^g and D ^ K^Iv^. If, for 1 < y < 1.5, f^'^^''Ky)/f ^^^'Ky) is 
almost independent of y then a precise value for \Vub\ can be extracted from the rates for 
B pii'e and B K*^!? integrated over this region in y (and /'^■^^^^(1)//'''^^^*''(1)). 
In a simple pole model this ratio of D decay form-factors is almost independent of y. We 
found that the matrix elements of the four-quark operators in the effective Hamiltonian for 
B K* II induce about a 20% uncertainty for the B —>■ K* il decay rate from cc threshold 
effects in the region 1 < y < 1.5. 

At the present time the rare decays B K* v v and B K* i i have not been observed, 
and there is no information on the individual form-factors for D —>■ plv^. Because of 
this, it is difficult to give a prognosis for the ultimate utility of the ideas presented here. 
However, even in the absence of the complete set of information needed for a high precision 
determination of \Vub\i our results may be useful. CLEO has reported a yield of about 1000 
B ^ p£i?£ events, corresponding to the branching ratio Br(i?° ivi) ~ (2 — 3) x 10~^ 

If heavy quark symmetry and SU{?>) are employed to get l/^'^^'^-'P (1 + 5^^^p'>) from 
the measured D — > K*!^^ form-factors, then eq. can be used to extract \Vub\ from the 
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large region of the Dalitz plot for the exclusive decay B — ^ pii^e. We predict, with this 
technique, a branching ratio of 5.2 {Vubl"^ for ii^e in the region 1 < y < 1.5. Lattice 



Monte Carlo simulations (and constituent quark model calculations ||21|) suggest that 
the violations of heavy quark symmetry and SU (3) symmetry that give corrections to the 
relation between f^^~*p^ and f(^~*^*) are not anomalously large. This method will give a 
value for \Vub\ that is on a more sound theoretical footing than that which results from a 
comparison of the endpoint of the electron spectrum of inclusive semileptonic B decay with 
phenomenological models. 

If experimental data on i? ^ K* u i? is available before a detailed study of semileptonic 
form-factors for D pii^i is performed, then using eq. (^) an extraction of |y('B^^*)|2 _j_ 
^(B^K*)^ should be possible. This gives a prediction for (l+S^^^'^^) with corrections 

of order m^, but no order l/rric corrections since heavy quark symmetry is not used. In this 
case there is no reason to restrict our analysis to the region of phase-space 1 < y < 1.5. 
Lattice QCD results suggest that the influence of SU{3) violation on the form-factors is 
small, and hence the value of \Vub\ that can be extracted in this way will be fairly precise. 
A sizable uncertainty in the theoretical prediction for the B K* i i decay rate arises 
from the charmonium resonance region. Without a better understanding of this, it will not 
be possible to extract \f(^~^^*)\'^ (1 + from this decay mode with high accuracy. 

Nonetheless, an extraction of (^i _|_ ^i^^K*)-^ from this mode may provide a useful 

determination of l/^'^^^-'P (1 + 5^^^^^) (and hence \Vub\) with uncertainties now from both 
SU (3) violation and from the contribution of the four-quark operators to the B — * K* £ i 
rate. 

Some improvements on the analysis in this paper are possible. Combining chiral per- 
turbation theory for mesons containing a heavy quark with heavy vector-meson chiral per- 



turbation theory allows a computation of the order m^lnms SU{3) violation in / p2[ . 
Unfortunately such an analysis cannot give a definitive result on the size of the SU{3) vi- 
olations because of unknown order counterterms. In this paper we have neglected the 
lepton masses. It is possible to include the corrections that arise from the non-zero value of 



the muon mass, although these are quite smalL 

A similar analysis to that performed in this paper can be done for the decays B ^ tt i ui, 
B^Kii, B^KuUjD—)- nlui and D Kiui. However, in these decays there are 
complications because very near zero recoil "pole contributions" I^S] spoil the simple scaling 
of the form-factors with the heavy quark mass. 
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